AO-A046  ^4b 


ONCLASSIFILD 


NAVAL  UNDERWATER  SYSTEMS  CENTER  NEW  LONDON  CONN  NEW  — ETC  F/6  20/14 
numerical  results  for  modified  IMAGE  THEORY  QUASI-STATIC  RANGE  — ETC(U) 
DEC  77  P R bannister*  R L DUBE 

NUSC-TR-5775  NL 


NUSC  Technical  Report  5775 

^Numerical  Results  for 
Modified  Image  Theory 
duasi  Static  Range 
Subsurf  aceto-Subsurf  ace 
and  Subsurface  to-Air 
Propagation  Equations 

P«t«r  R.  Bannister 
Ran#'  L.  Duba 
Sub  marina  Eiectromagnatic 
Systems  Dapartmant 


D n c 


'*  7 3 i97fi 

I i 

ijLbiii;LbO  l/1^ 


7 Daeambar  1977 


NUSC 


NAVAL  UNOERWATIR  SYSTEMS  CENTER 
Nawport.Rhods  ■•Und  • N«ww  London.Conn*etiout 


Approved  for  public  release;  distribution  unlimited 


PREFACE 

This  report  was  prepared  under  NUSC 
Project  No.  A-532-24,  P.  R.  Bannister, 
Principal  Investigator;  Naval  Industrial 
Funding  was  from  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center,  W,  J. 
Andahazy,  Project  Director.  The  sponsor 
was  the  Naval  Sea  Systems  Command,  Project 
Program  Manager,  W.  L.  Welsh,  SEA  03424. 

The  Technical  Reviewer  for  this  report 
was  J.  J.  Tennyson  (Code  3403). 


REVIEWED  AND  APPROVED:  7 December  1977 


John  Merrill 

Head:  Submarine  Electromagnetic 
* Systems  Department 


The  authors  of  this  report  are  located  at  the  New  London 
Laboratory,  Naval  Underwater  Systems  Center, 

New  London,  Connecticut  06320. 


REPORT  DOCUMENTATION  PAGE 

READ  DtSTRUCTIONS 

BEFORE  COMPLETWC  FORM 

1.  MlfoAT  NUmIIa  / 2.  OOVT  ACCKUION  NO. 

TR  5775 

. 7=" 

i.  RECIRIENT'S  CATALOG  NUMBER 

^ NUMERICAL,JIESULTS  FOR,  MODIFIED.  IMAGE, THEORY 
^'^lUAS  I - STATIC., RANGE . SUBSURFACE- f6:;SUBSURFACE 
AND  SUBSURFACE-TO-AIR  .PROPAGATION^ EQUATIONS r 

\ ! 

' Pete^R. /Bannister 

Ren^L.^Dube  ^ 

■■■ 

0.  NCRrONMINO  ONOANIZATION  NAME  ANO  AODNISt 

Naval  Underwater  Systems  Center 

New  London  Laboratory 

New  London.  CT  06.'?20 

10.  RRQGRAM  ELEMENT.  RROJECT.  TASK 
AREA  A WORK  UNIT  NUMBERS 

A-532-24 

It.  CONTAOULINO  OFTlCt  HANE  AMO  AODREt*  /^i  ! 

Naval  Sea  Systems  Command  ^ U j 

SEA-03424 

Washington,  DC  20362 

IS.  NUMBER  OF  RAGES 

30 

U.  monitoring  agency  name  O AOOREESfll  Mllaianl  Iram  CanHalllnt  OlHea) 

If.  fECUmTY  CLASS.  (9t  Ml*  raporfj 

UNCLASSIFIED 

Rif.i  1 1 1 1 1 1 11 1 mm 

!«.  OitTRIRUTION  STATEMENT  (at  IMa  Kaparl) 

Approved  for  public  release;  distrj.bution  unlimited. 


IT.  CMITNiaUTIOM  tTATeiUMT  (•!  »*•  mifd  In  ViMk  JO.  II  OIWnMnI  htm  lt«onM> 


10.  SUP^LIMINTARY  MOTCf 


10.  KBV  WOAOt  fCnnilnun  an  rarataa  ml4a  II  naeaaaaar  an^  HanMIy  Or  WneO  niaiOrr; 


to. 


Dipole  Antennas 
Modified  Image  Theory 
iasi-Static  Range 
Niimerical  Integration  Results 


Subsurface- to-Air  Propagation 
Subsurface- to-Subsur face  Propagation 
Earth -Image  Theory 


(Canllnua  an  taaataa  alda  II  naaaaaair  amd  I4anlllt  Or  OInoO  niwOwj 

Numerical  results  are  presented  for  the  previously  derived  horizontal 
electric  dipole  modified  image-theory  quasi-static  range  field  component 
approximations.  Both  subsurface- to- subsurface  and  subsurface-to-air  propa- 
gation cases  are  considered.  It  is  shown  that  the  modified  image  theory 
approximations  are  in  good  agreement  with  the  more  exact  numerical  integration 
results. 


K 


TABLE  OF  CONTENTS 


TR  5775 


LIST  OF  ILLUSTRATIONS 
INTRODUCTION  .... 


HORIZONTAL  ELECTRIC  DIPOLE  (HED)  FIELD  COMPONENT 
EXPRESSIONS  


NUMERICAL  RESULTS 
CONCLUSIONS  . . . 
REFERENCES  . . . 


JiW'S;.  1 ' 

hV. 

I 

( 

i 

j ii‘.;rsn 

j li’:l 

- 

D O r 


^ iHi 

Jr-  IA:-!  13  1978  I 

D 


TR  5775 


; 

1 


i 

[f 


I 


LIST  OF  ILLUSTRATIONS 


Figure 

1 Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  Magnetic  Fields  in 
Air  Produced  by  a Buried  HED  (H  = -1.0,  Z = 1.0)  . . 

2 Comparison  of  Composite  Modified  Image  Theory  and 

Numerical  Integration  Results  for  the  Magnetic 

Fields  in  Air  Produced  by  a Buried  HED 

CH  = -1.0,  Z = 1.0) 

3 Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  Ep  Component  in  Air 
Produced  by  a Buried  HED  (H=-2.0,  Z=1.0).  . . . 

4 Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  El  Component  in  Air 
Produced  by  a Buried  HED  (H  = -2.0,  Z = 1 . 0)  . . . . 

5 Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  EJ  Component  in  Air 
Produced  by  a Buried  HED  (H  = -2.0,  Z = 1 . 0)  . . . . 

6 Comparison  of  Modified  Image  Theor>'  and  Numerical 

Integration  Results  for  the  Ep  Component  in  Air 
Produced  by  a Buried  HED  (H=-1.0,  Z=2.0)  . . . . 

7 Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  El  Component  in  Air 
Produced  by  a Buried  HED  (H  = -1.0,  Z = 2.0)  . . . . 

8 Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  Ej  Component  in  Air 
Produced  by  a Buried  HED  (H  = -1.0,  Z = 2.0)  . . . . 

9 Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  Hp  Component  in  Air 
Produced  by  a Buried  HED  (H  = -1.0,  Z = 2.0)  . . . . 

10  Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  HI  Component  in  Air 
Produced  by  a Buried  HED  (H=-1.0,  Z=2.0)  , . . . 

11  Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  HJ  Component  in  Air 
Produced  by  a Buried  HED  (H  = -1.0,  Z = 2.0)  . . . . 

12  Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  Hp  Component  in  Air 
Produced  by  a Buried  HED  (H  = -2 . 0,  Z =■  1 . 0)  . . . . 

13  Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  H^  Component  in  Air 
Produced  by  a Buried  HED  (H  = -2 . 0,  Z = 1 . 0)  . . . . 


Page 


9 


10 

11 

11 


12 


12 

13 

13 

14 

14 

15 

15 

16 


I'! 

ji 


j 

i 


I 


ii 


Figure 


TR  5775 


LIST  OF  ILLUSTRATIONS  (Cont'd) 


Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  H$  Component  in  Air 
Produced  by  a Buried  HED  (H  = -2.0,  Z = 1.0)  . . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  Ep  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -1.0)  . . . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  El  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -1.0)  . . . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  Hp  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -1.0)  . . . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  H^  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -1.0)  . . . . 

Comparison  of  Modified  Image  Theor>'  and  Numerical 

Integration  Results  for  the  H£  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -1.0)  . . . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  Hp  Component  in  Earth 
Produced  by  a Buried  HED  (H  = -1.0,  Z = 0)  . . . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  H|  Component  in  Earth 
Produced  by  a Buried  HED  (H  = -1.0,  Z = 0)  . . . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  Hp  Component  in  Earth 
Produced  by  a Buried  HED  (H  = -1.0,  Z - -1.0)  . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  H^  Component  in  Earth 
Produced  by  a Buried  HED  (H  = -1 . 0,  Z = - 1 .0)  . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  Ep  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -2.0)  . . . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  E^  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -2.0)  . . . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  Hp  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -2.0)  . . . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  H^  Component  in  Earth 
Produced  by  a Buried  HED  fH  = 0,  Z = -2.0)  . . . . 

Comparison  of  Modified  Image  Theory  and  Numerical 

Integration  Results  for  the  Hi  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -2.0)  . . . . 


iii/iv 
Reverse  Blank 


TR  5775 


NUMERICAL  RESULTS  FOR  MODIFIED  IMAGE  THEORY 
QUASI -STATIC  RANGE  SUBSURF ACE-TO-SUBSURFACE 
AND  SUBSURFACE-TO-AIR  PROPAGATION  EQUATIONS 


INTRODUCTION 


During  the  past  few  years,  considerable  interest  has  been  shown  in 
determining  the  quasi-static  field  components  of  antennas  located  above 
or  buried  beneath  the  earth’s  surface.  The  quasi-static  range  is  de- 
fined as  the  range  where  the  measurement  distance  is  much  less  than  a 
free-space  wavelength.  Quasi-static  range  results  are  useful  for  sub- 
marine radio  communication  and  detection  as  well  as  for  the  buried-miner 
problem.  They  are  also  helpful  to  geophysicists  engaged  in  determining 
the  electrical  properties  of  the  earth. 

Some  work  has  been  done^"®  on  determining  the  quasi-static  fields 
from  various  subsurface  sources  when  the  measurement  distance  R = 

(p2  + z2^1/2  comparable  to  the  earth  skin  depth  6.  However,  the 
resulting  field  strength  expressions  are  very  complicated  because  they 
involve  products  of  modified  Bessel  functions  of  different  argument. 
Computing  the  field  component  expressions  is  lengthy  and  difficult,  but 
some  numerical  results  have  been  obtained.  One  method  of  obtaining 
these  results  has  been  discussed  by  Atzinger,  Pensa,  and  Pigott.®  Nu- 
merical integration  techniques  have  also  been  employed.^'® 

Recently,®  by  utilizing  techniques  of  finitely  conducting  earth 
image  theory,  the  authors  have  derived  approximate  expressions  for  the 
general  quasi- static  range  electromagnetic  fields  produced  by  various 
subsurface  antennas.  It  is  the  purpose  of  this  report  to  provide  some 
numerical  calculations  of  these  approximate  field  component  expressions. 
The  image-theory  expressions  are  also  compared  with  the  more  exact 
numerical  integration  results. 

Physically,  the  essence  of  the  finitely  conducting  earth-image 
theory  technique  is  to  replace  the  finitely  conducting  earth  by  a 
perfectly  conducting  earth  located  at  the  (complex)  depth  d/2,  where 
d = 2/y  = 6(l-i).  Analytically,  this  corresponds  to  replacing  the 
algebraic  reflection  coefficient  (u-A)/ (u+A)  in  the  exact  integral 
equations  by  exp(-Ad),  where  A is  the  variable  of  integration. 

For  antennas  located  at  or  above  the  earth's  surface,  the  general 
image  theory  approximation  is  valid  throughout  the  quasi-static 
range. However,  if  one  or  both  antennas  are  buried,  the  pre- 
viously derived^®  surface-to-surface  and  surface-to-air  image-theory 
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results,  multiplied  by  the  exponential  attenuation-with- depth  factor 
exp[Y(z  + h) ] or  expCyh),  generally  will  be  valid  only  for  R > 3|z  + h|, 
where  h and  z,  respectively,  are  the  depth  of  the  source  and  receiving 
antennas. 

Therefore,  we  modified  further  the  exact  subsurface- to- subsurface 
and  subsurface- to-air  integral  equations  and  obtained  fairly  simple 
field- strength  equations^  for  the  general  quasi-static  range,  which 
reduce  to  previously  derived  results  when  R » 6 and  when  R <<  6 and 
|z  + h|  « 6. 

For  the  purposes  of  this  report,  we  will  consider  a horizontal 
electric  dipole  (HED)  source  that  is  oriented  in  the  x-direction  an^ 
located  at  depth  h(h  ^ 0)  with  respect  to  a cylindrical  coordinate  sys- 
tem (p,  (t>,  z)  . The  dipole  is  assumed  to  be  of  infinitesimal  length  Z 
and  carry  a constant  current  I.  The  earth  occupies  the  lower  half-space 
(z  < 0)  and  the  air  occupies  the  upper  half-space  (z  > 0).  Displace- 
ment currents  are  neglected  in  both  the  ground  and  the  air.  The  mag- 
netic permeability  of  the  earth  is  assumed  to  equal  the  permeability  of 
free  space  Po-  Meter-kilogram- second  (MKS)  units  are  employed  and  a 
suppressed  time  factor  of  exp(ia)t)  is  assumed. 


HORIZONTAL  ELECTRIC  DIPOLE  (HED)  FIELD 
COMPONENT  EXPRESSIONS 


The  appropriate  modified  image-theory  equations  for  HED  subsurface- 
to-air-propagation  (h  < 0,  z > 0)  are  given  on  page  6 of  the  previous 
report  ^ and  are  as  follows: 
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^ z 4ir 

j where 

K^  = p2  + (2  . bh)2  and  K|  = p2  + (d  + z - bh)^  , 

' a = 0 and  b = 1 for  R <<  6 and  |h|  <<  6, 

a = 0.4  and  b =0.96  for  R/6  less  than  approximately  1, 

< a = 0.96  and  b = 0.4  for  R/6  between  approximately  1 and  10,  and 

a = 1.0  and  b = 0 for  R > |3h|. 

i 

I The  appropriate  modified  image  theory  equations  for  HED  subsurface- 

i to-subsurface  propagation  (h  < 0,  z < 0) , given  on  pages  13  and  14  of 

the  previous  report,^  are 
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where 


K3  = + [b(z  + h]^  and  + [d  - b(z  + h)]^  , 

Rg  = P^  + (z  - h)^  and  r|  = P^  + (z  + h)^, 
a = 0 and  b = 1 for  Ri/6  <<  1, 

a = 0.4  and  b = 0.96  for  Ri/6  less  than  approximately  1, 
a = 0.96  and  b = 0.4  for  Ri/6  between  approximately  1 and  10,  and 
a = 1 and  b = 0 for  p > 3|z  + h|. 


NUMERICAL  RESULTS 


Note  that  we  have  defined  loosely  the  variables  a and  b.  That  is, 
we  let  a = 0.96  and  b = 0.4  for  Ri/6  between  approximately  1 and  10,  and 
a = 0.4  and  b = 0.96  for  Ri/6  less  than  approximately  1.  The  specific 
crossover  point  for  each  field-strength  component  will  depend  not  only 
on  Ri/6  but  also  on  | z + h|/6. 

Because  the  resulting  field-strength  formulas  can  be  calculated 
easily  on  a dosk-top  calculator,  the  field  strengths  can  be  determined 
by  using  both  values  of  a and  b,  thus  numerically  determining  the  cross- 
over point.  Alternatively,  we  could  set  the  two  expressions  (involving 
different  values  of  a and  b)  equal  and  solve  for  the  crossover  point. 

Comparisons  of  modified  image  theor>'  and  numerical  integration 
results  for  the  electric  and  magnetic  fields  in  air  (or  earth)  pro- 
duced by  a buried  HED  are  presented  in  figures  1 through  28,  which  fol- 
low page  8.  The  normalized  amplitude  of  each  component  (E'  or  H')  is 
plotted  versus  p/6  for  various  values  of  H(=  h/6)  and  Z(=  z/6) , where 


ir  = 
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4Tr6^H 
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(13) 


The  normalization  factor  (0  dB)  is  the  numerical  integration  value 
of  each  component  at  p/6  = 0.1. 

Presented  in  figures  1 and  2 are  comparisons  of  the  modified  image 
theory  and  numerical  integration  results  for  the  magnetic  fields  at  a 
height  of  one  skin  depth  (Z  = 1.0)  produced  by  an  HED  buried  at  a depth 
of  one  skin  depth  (H  = -1.0).  From  figure  1,  we  see  that  the  crossover 
point  for  this  example  is  p/6  = 1.5  (R/6  = 1.8).  That  is,  if  0.1  < 
p/5  < 1.5,  then  a = 0.4  and  b = 0.96,  whereas  if  1.5  < p/6  < 10,  then 
a = 0.96  and  b = 0.4.  A comparison  (figure  2)  of  the  composite  modi- 
fied image-theory  results  with  the  numerical  integration  results  shows 
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that  the  modified  image-theory  calculations  are  within  approximately 
1 dB  of  the  numerical  integration  results  over  the  complete  range  of 
p/6  (0.1  to  10).  In  fact,  the  modified  image  theory  even  predicts  the 
Hp  component  amplitude  dip  at  the  right  place  Cp/6  - 1.5). 

Comparisons  of  the  modified  image-theory  and  numerical  integration 
results  for  the  electric  fields  at  a height  of  one  skin  depth  produced 
by  an  HED  buried  at  two  skin  depths  (H  = -2.0)  are  presented  in  figures 

3 through  5.  From  figure  3,  we  see  that  the  a = 0.96,  b = 0.4  curve  is 

in  better  agreement  with  the  numerical  integration  result  throughout 
the  range  of  p/6  plotted.  For  values  of  p/6  > 3,  a much  better  fit  to 
the  numerical  integration  re.sult  can  be  obtained  by  letting  a = 1 and 
b = 0.  Alternatively,  one  could  employ  the  quasi-near  range  formulas 
given  in  table  3.13  of  Kraichman.^^ 

From  figure  4,  we  see  that  the  a = 0.96,  b = 0.4  curve  is  in  good 
agreement  with  the  numerical  integration  result  for  p/6  = 0.1  to  10. 
However,  for  the  E^  component  (figure  5),  the  crossover  point  is  at 
p/6  ' 1.5.  That  is,  if  0.1  < p/6  < 1.5,  then  a = 0.4  and  b = 0.96, 

whereas  if  1.5  < p/6  < 10,  then  a = 0.96  and  b = 0.4. 

Presented  in  figures  6 through  8 are  comparisons  of  the  modified 
image-theory  and  numerical  integration  results  for  the  electric  fields 
at  a height  of  two  skin  depths  produced  by  an  HED  buried  at  one  skin 
dejith.  From  these  figures  we  observe  that  the  a = 0.96,  b = 0.4  curve 
is  in  excellent  agreement  with  the  numerical  integration  result  through- 
out the  range  of  p/6  considered  (0.1  to  10). 

It  should  be  noted  that  because  of  reciprocity,  the  normalized 
electric  field  plots  presented  in  figures  3 through  8 are  also  valid 
for  air-to-subsurface  propagation  if  Z and  H are  interchanged.  That  is, 
figures  3 and  4 also  represent  the  horizontal  electric  field  at  a depth 
of  two  skin  depth  produced  by  a horizontal  electric  dipole  located  at  a 
height  of  one  skin  depth,  while  figure  5 represents  the  horizontal 
electric  field  at  a depth  of  two  skin  depths  produced  by  a vertical 
electric  dipole  (VED)  source  at  a height  of  one  skin  depth.  For  fur- 
ther details,  see  reference  9. 

Comparisons  of  the  modified  image  theory  and  numerical  integration 
results  for  the  magnetic  fields  at  a height  of  two  skin  depths  produced 
by  an  HED  buried  at  one  skin  depth  are  presented  in  figures  9 through  11. 
From  figure  9,  we  see  that  the  crossover  point  for  the  Ifp  component  is 
at  p/6  ' 2.2.  That  is,  if  0.1  < p/6  < 2.2,  then  a = 0.4  and  b = 0.96, 
whereas  if  2.2  < p/6  < 10,  then  a = 0.96  and  b = 0.4.  For  the  H^  com- 
ponent (figure  10),  the  a = 0.96,  b = 0.4  curve  is  in  excellent  agree- 
ment with  the  numerical  integration  result  throughout  the  range  of  p/6 
considered.  From  figure  11,  we  see  that  the  crossover  point  for  the  H^ 
component  is  at  p/6  - 1.5. 
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Presented  in  figures  12  through  14  are  comparisons  of  the  modified 
image  theory  and  numerical  integration  results  for  the  magnetic  fields 
at  a height  of  one  skin  depth  produced  by  a buried  HED  at  a depth  of 
two  skin  depths.  From  these  figures,  we  see  that  the  crossover  point 
is  p/6  - 2 for  the  H'  component  (figure  12),  p/6  > 0.8  for  the  HI  com- 
ponent (figure  13),  and  p/6  - 1.7  for  the  component  (figure  14). 

It  should  be  noted  that  because  of  reciprocity,  the  normalized 
magnetic  field  plots  presented  in  figures  9 through  14  are  also  valid 
for  air-to- subsurface  propagation  if  Z and  H are  interchanged.  That  is, 
figures  9 and  10  also  represent  the  horizontal  electric  field  at  a 
depth  of  one  skin  depth  produced  by  a horizontal  magnetic  dipole  (HMD) 
source  located  at  a height  of  two  skin  depths  while  figure  11  repre- 
sents the  horizontal  electric  field  at  a depth  of  one  skin  depth  pro- 
duced by  a vertical  magnetic  dipole  (VMD)  source  located  at  a height  of 
two  skin  depths.  For  further  details,  see  reference  9. 

Comparisons  of  the  modified  image  theory  and  numerical  integration 
results  for  the  electric  and  magnetic  fields  at  a depth  of  one  skin 
depth  produced  by  an  HED  located  at  the  earth's  surface  (H=0)  are  given 
in  figures  15  through  19.  From  these  figures,  we  see  that  the  cross- 
over point  is  p/6  - 1.5  for  the  Ep  component  (figure  IS),  p/6  > 1.0  for 
the  E|  component  (figure  16),  p/6  - 1.5  for  the  Up  component  (figure 
17),  p/6  - 1.0  for  the  H^  component  (figure  18),  and  p/6  - 1.0  for  the 
component  (figure  19).  It  should  be  noted  that  for  p/6  > 3,  an  even 
better  fit  to  the  numerical  integration  results  can  be  obtained  by 
setting  a = 1.0  and  b = 0.  Alternatively,  one  could  employ  the  quasi- 
near  range  formulas  given  in  table  3.16  of  Kraichman.^^ 

Referring  to  equations  (7)  through  (11),  we  see  that,  for  the  Ep, 
E^,  and  H^  components,  identical  results  will  be  obtained  for  the  H=0, 
Z=-l  and  H=-l,  Z=0  cases.  However,  for  the  Hp  and  H^  components,  this 
will  not  be  the  case.  This  fact  is  further  illustrated  in  figures  20 
and  21,  which  present  comparisons  of  the  modified  image-theory  and 
numerical  integration  results  for  the  horizontal  magnetic  fields  at  the 
surface  (Z=0)  produced  by  a HED  buried  at  a depth  of  one  skin  depth 
(H=-l).  From  figure  20,  we  see  that  the  a = 0.96,  b = 0.4  Hp  curve  is 
in  fair  agreement  with  the  numerical  integration  result  for  0.1  < 
p/6  < 3 and  in  very  good  agreement  for  p/6  > 3.  Referring  to  figure 
21,  we  see  that  the  crossover  point  for  the  H,j,  component  is  at  p/6  - 
1.5.  That  is,  if  0.1  < p/6  < 1.5,  then  a = 0.4  and  t = 0.96,  whereas 
if  1.5  < p/6  < 10,  then  a = 0.96  and  b = 0.4. 

Presented  in  figures  22  and  23  are  comparisons  of  the  modified 
image  theory  and  numerical  integration  results  for  the  horizontal  mag- 
netic fields  at  a depth  of  one  skin  depth  produced  by  a HED  buried  at 
one  skin  depth.  From  these  figures  we  see  that  the  crossover  point  for 
both  components  occurs  at  p/6  1.0.  It  should  be  noted,  however,  that 
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for  the  component  (figure  23),  the  a = 0.96,  b = 0.4  curve  is  in 
good  agreement  with  the  numerical  integration  result  throughout  the 
range  of  p/6  considered. 

Comparisons  of  the  modified  image  theor>'  and  numerical  integration 
results  for  the  electric  and  magnetic  fields  at  a depth  of  two  skin 
depths,  produced  by  an  HED  located  at  the  earth's  surface,  are  presented 
in  figures  24  through  28.  From  these  figures  we  see  that  the  crossover 
point  is  at  p/6  - 1 to  1.5  for  all  five  components.  The  agreement  be- 
tween the  modified  image  theory  and  numerical  integration  results  is 
very  good  for  the  Hp  and  components  (figures  26  and  27)  and  fair  for 
the  Ep,  E^,  and  components  (figures  24,  25,  and  28).  The  modified 
image  theory  fails  to  predict  the  Ep  amplitude  dip  at  p/6  - 2.5  (figure 
24),  although  it  did  predict  the  Ep  amplitude  dip  for  the  H=0,  Z=-l 
case  (figure  15).  For  both  the  E^  and  components  (figures  25  and 
28),  the  a=0.4,  b=0.96  curve  is  in  good  agreement  with  the  numerical 
integration  result  from  0.1  < p/6  < 3.  For  p/6  > 3,  better  results  can 
be  obtained  by  setting  a=1.0  and  b=0. 

Because  of  reciprocity,  the  normalized  electric  and  magnetic  field 
plots  presented  in  figures  15  through  28  are  also  valid  for  subsurface 
(or  surface)  to  subsurface  (or  surface)  propagation  if  Z and  H are 
interchanged.  That  is,  figures  15  and  16  also  represent  the  horizontal 
electric  fields  at  the  surface  produced  by  an  HED  buried  at  one  skin 
depth;  figures  17  and  18  also  represent  the  horizontal  electric  fields 
at  the  surface  produced  by  an  HMD  buried  at  one  skin  depth;  and  figure 
19  also  depicts  the  horizontal  electric  field  at  the  surface  produced 
by  a buried  VM)  located  at  H=-l.  Reference  9 contains  further  details. 


CONCLUSIONS 

Numerical  results  have  been  presented  for  the  previously  derived^ 
HED  modified  image-theory  quasi-static  range  subsurface- to- subsurface 
and  subsurface- to-air  propagation  approximate  field  component  expres- 
sions. It  has  been  shown  that,  for  the  most  part,  the  modified  image- 
theory  approximations  are  in  good  agreement  with  the  more  exact  numerical 
integration  results. 

These  modified  image  theory  results  will  be  useful  for  submarine 
radio  communication  and  detection  as  well  as  for  the  buried  m.iner  prob- 
lem. They  also  may  be  helpful  to  geophysicists  engaged  in  determining 
the  electrical  properties  of  the  earth. 
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Figure  2.  Comparison  of  Composite  Modified  Image  Theory  and 
Numerical  Integration  Results  for  the  Magnetic  Fields  in  Air 
Produced  by  a Buried  HED  (H  = -1.0,  Z = 1.0) 
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Figure  3.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  Ep  Component  in  Air  Produced  by  a 
Buried  HED  (H  = -2.0,  Z = 1.0) 


Figure  4.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  E^  Component  in  Air 
Produced  by  a Buried  HED  (H  = -2.0,  Z = 1.0) 
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Figure 


5.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  E2  Component  in  Air 
Produced  by  a Buried  HED  (H  = -2.0,  Z = 1.0) 


6.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  E'  Component  in  Air 
Produced  by  a Buried  HED  (H  = -1.0,  Z = 2.0) 
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Figure  7.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  El  Component  in  Air 
Produced  by  a Buried  HED  (H  = -1.0,  Z = 2.0) 
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Figure  8.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  Ej  Component  in  Air 
Produced  by  a Buried  HED  (H=  -1.0,  Z = 2.0) 
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Figure  9.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  Hp  Component  in  Air 
Produced  by  a Buried  HED  (H  = -1.0,  Z = 2.0) 


Figure  10.  Comparison  of  Modified  Image  Theor>-  and  Numerical 
Integration  Results  for  the  Component  in  Air 
Produced  by  a Buried  HED  (H  = -1.0,  Z = 2.0) 
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I Figure  13,  Comparison  of  Modified  Image  Theory  and  Numerical 

J Integration  Results  for  the  Component  in  Air 

Produced  by  a Buried  HED  (H  = -2.0,  Z = 1.0) 


Figure  14.  Comparison  of  Modified  Image  Theor}'  and  Numerical 
Integration  Results  for  the  Component  in  Air 
Produced  by  a Buried  HED  (H  = -2.0,  7 = 1.0) 
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Figure  15.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  £„  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -1.0) 
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Figure  16.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  El  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -1.0) 
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Figure  17.  Comparison  of  Modified  Image  Theory  and  Num.erical 
Integration  Results  for  the  H'  Component  in  Farth 
Produced  by  a Buried  HED  (H  = 0,  Z = -1.0) 


Figure  18.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -1.0) 
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21.  Comparison  of  Modified  Image  Theory  and  Numerical 
[ntegration  Results  for  the  H$  Component  in  Earth 
Produced  by  a Buried  HED  (H  = -1.0,  Z = 0) 


Figure  22.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  Up  Component  in  Earth 
Produced  by  a Buried  HED  (H  = -1.0,  Z = -1.0) 
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Figure  23.  Comparison  of  Modified  Image  Theor>-  and  Numerical 
Integration  Results  for  the  Component  in  Earth 
Produced  by  a Buried  HED  (H  = -1.0,  Z = -1.0) 


Figure  24.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  Ep  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -2.0) 
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25.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -2.0) 
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Figure  26.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  Hp  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -2.0) 
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Figure  27.  Comparison  of  Modified  Image  Theory  and  Numerical 
Integration  Results  for  the  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -2.0) 


Figure  28.  Comparison  of  Modified  Image  Theory-  and  Numerical 
Integration  Results  for  the  H2  Component  in  Earth 
Produced  by  a Buried  HED  (H  = 0,  Z = -2.0) 
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